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ABSTRACT 

We present an analysis of the X-ray emission of a large sample of z 3 Lyman break 
galaxies (LBGs), based on Chandra/ ACIS observations of several LBG sm-vey fields. A 
total of twenty-four LBGs are directly detected in the X-ray, approximately doubling 
the number of known detections. Thirteen of the LBGs have optical spectroscopic 
signatures of AGN activity, but almost all the other X-ray detections are also likely to 
host an accreting black hole based on their X-ray properties. The AGN exhibit a wide 
range in X-ray luminosity, from weak Seyferts to bright QSOs. Optical spectroscopy 
identified approximately 1/3 of the X-ray detected sources as broad hue QSOs, 1/3 
as narrow line AGN and 1/3 as normal star forming LBGs. The fraction of X-ray 
detected LBGs is 3 per cent, much lower than has been found for submm selected 
galaxies. Two galaxies have X-ray luminosities, spectra and fx/fopt values that are 
consistent with emission from star formation processes and are identified as candidate 
X-ray bright, pure starburst galaxies at z ~ 3. If powered solely by star formation the 
sources would have SFRs of 300-500 M0 yr~^. X-ray spectral analysis of the LBGs 
shows a mean photon index of F = 1.96, similar to local AGN. There is evidence 
for absorption in at least 40 per cent of the objects. Significantly more absorption is 
evident in the narrow line AGN, which is consistent with AGN unification schemes. 
After correction for absorption the narrow and broad line objects show the same av- 
erage luminosity. X-ray detected LBGs spectroscopically classified as normal galaxies, 
however, are less luminous in both soft and hard X-ray bands indicating that the host 
galaxy is outshining any optical AGN signature. Turning to the X-ray emission from 
LBGs without direct detections, stacking the X-ray flux in the two deepest Chandra 
fields under consideration (the HDF-N and GWS) produced significant detections in 
each, although the GWS result was marginal. The detection in the HDF-N gives an 
X-ray derived SFR of 42.4 ± T.SMq yr~^ per LBG and, by comparing with the UV 
SFR, the implied UV extinction correction is 4.1 ±0.8. The LBG sample was split into 
three bins based on UV magnitude to examine the correlation between UV and X-ray 
emission: for the limited statistics available there was no evidence of any correlation. 

Key words: galaxies: active - galaxies: starburst - galaxies: high-redshift - X-rays: 
galaxies 



1 INTRODUCTION 

Observations of distant galaxies are an essential compo- 
nent in our understanding of the development of the Uni- 
verse. Most high-redshift galaxies have b een identified in the 
UV via the Lyman break technique (e.g. JSteidel et alJll996t 
iLowenthal et al.lll997l: ISteidel et alj |2003h. There is there- 
fore considerable interest in constraining properties of LBGs 
such as the star- formation rate (SFR), dust content, stellar 
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mass and metallicity. Another effective way of identifying 
high-z galaxies is by observations at sub-mm wavelengths 
with SCUB A on the James Clerk Maxwell Telescope (e.g. 
Ilvison et a l. 2002; Smail ct al. 2002; Chapman et al. 200$. 
The relation between LBGs and SCUBA galaxies is still not 
entirely clear but studies have shown that up to 50 per cent 
of SCU BA galaxies have s imilar rest-frame U V colours as 
LBGs jReddv et al.ll2006l: ISteidel et alJl2004ll . LBGs and 
SCUBA galax;ies may simply form a continuous distribution 
in SFR and dust content, with S CUBA galaxies occupying 
the upper end of the distribution llReddv et alJl2006l . 12003 : 
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lAdelbereer fc Steidelllioooll . Indeed a recent study of IR lu- 
minous LBGs showed that while LBGs have a diverse range 
of IR properties, at least some LBGs have the same IR prop- 
erties as SCUBA galaxies and are li kely to be the progen itors 
of today's massive giant ellipticals llHuang et al.ll2005l) . 

The proposed starburst-AGN connection that arises 
out of merger-driven galaxy formation scenarios (e.g. 
Sanders fc Mi rabel 1996; H opkins et alJl200a) . and which is 
a natural explanation for the observed Mbh — u relation 
between the mass of dormant black holes and galaxy bulges 
in local galaxies f e.g.lMag orrian et alJll993 : lGebhardt et alJ 
[2000 ; Fcrrarcso fc Merritd I2OO0I) . also leads to an interest 
in the AGN content of high-z star forming galaxies. The 
incidence of AGN in SC UBA sources has attracted con- 
sider able attention (e.g. iFab ian et aL J2000|; iBarger et alJ 
I2OOII) . Most recentlv. lAlexander et alj ^200^ have found an 
extremely high rate of X-ray detection of radio-identified 
SCUBA galaxies, inferring that at least 40 per cent, and 
perhaps as many as 75 per cent of submm-selected galax- 
ies contain an AGN. Much less work has been done on the 
properties of AGN in the LBG population. To date most 
work including X-ray data has focussed on SF Rs and extinc- 
tion corrections of LBGs without AGN (e.g . iBran^^^^ 
l200ll:lNandra et al.ll2002l:lLehmer et al.ll2005n . lSteidel et~ 
l|2nn2|^ and lHunt et alJ lj2004) analysed rest-frame UV spec- 
tra of LBGs to identify AGN based on emission lines and to 
calculate the faint end of the AGN luminosity function at 
high-2. They found that 3 per cent of LBGs harbour AGN, 
which contribute 8 p er cent to the integrated 1500 A UV 
luminosity at 2 = 3. iNandra. Laird, fc Steidel (|2QQ^ used 
LBG selection in two fields (the HDF-N and GWS) as a 
method of efficiently identifying the redshift of AGN in X- 
ray surveys. The results were used to estimate the space 
density of moderate luminosity AGN at 2: = 3, finding them 
to be ten times more common than high luminosity QSOs. 

In this paper the X-ray emiss ion of z ^ 3 LBGs i n six 
of the survey fields presented by ISteidel et alj (l2003l) are 
analyzed. This work extends the analysis of LBGs in the 
1 Ms CDF-N by N andra et al. (2002; hereafter N02) and 
iNandra. Laird, fc S teidel (2005) by using more LBG fields 
and including X-ray spectral analysis of the LBGs. Follow- 
ing on from the s tacking analys i s of L BGs bv lBrandt et alJ 
J200ll) . N02 and iLehmer et al.l (1200511 . the X-ray emission 
from non-AGN dominated LBGs in the two deepest fields is 
also analyzed. 

Throughout, a standard, fiat ACDM cosmology with 
JIa ~ 0.7 and Ho — 70 km s~^ Mpc^^ is assumed. 



2 DATA AND SOURCE SAMPLE 
2.1 Optical LBG data 

Our s ample of 2~3 LBGs is culled from the ISteidel et alJ 
l|2003() survey fields that have archived Chandra ACIS imag- 
ing data. These fields are the Hubble Deep Field-North 
(HDF-N), Groth-Westphal Strip (GWS), Q1422+2309, 
SSA22 (SSA22a and SSA22b), B20902-H34 and 3C 324. 

The areas of overlap between the LBG fields and the 
Chandra data, along with the number of LBGs in the over- 
lap areas are shown in Table [T] The LBG data for fields 
B20902-f 34 and 3C 324, for which the targets of both were 
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Figure 1. Survey area versus tlie minimum effective exposure in 
the soft band exposure maps for all of the LBG fields included 
in this work. The hatched and shaded sections show the survey 
area versus minimum effective exposure for the LBG regions in 
the GWS and HDF-N fields, respectively. 



radio galaxies, covers a small area and therefore contain rela- 
tively few LBG candidates. Q1422-I-2309, which was centred 
on a z = 3.620 gravitational l y lens ed QSO, has the deepest 
data of all the lSteidel et al.l ll2003li survey fields. 

In the HDF-N region, imaging for LBG selection and 
spectroscopy has recently been extended by the Steidel 
team from the original survey area to now cover the larger 
GOODS-N region. The catalogue of HDF-N LBGs used in 
this paper consists of all the LBGS from lsteidel et al.l (l20o3l 
plus the LBGs from the extended survey area. The latest 
imaging and photometry is of better quality than the orig- 
inal data and is therefore used for any LBG candidates in 
the original list that are also covered by the new data. 

Of the LBGs considered in this work 43 per cent have 
confirmed spectroscopic identifications as 2 > 2 galaxies. 
Over the whole LBG survey (containing more than 2300 
objects) the redshift distribution of candidates is 2: = 2.96 ± 
0.29 and the fraction of low-redshift interlopers is only 4 per 
cent. Given the tightly peaked redshift distribution and the 
low interloper fraction we are confident that the remainder 
of our sample without spectroscopic confirmation are indeed 
high-redshift star-forming galaxies. 



2.2 X-ray data and reduction 

The majority of the results in this paper come from the 
deep X-ray observations of the HDF-N and GWS, with the 
remaining fields acting as a supplementary data set. The 
Chandra ACIS-I observation of the HDF-N (the Chandra 
Deep Field-North, CDF-N), with a total exposure time of 
approximately 2 Ms, is the deepest X-ray observation taken 
to date. A full description of the data and point source anal- 
ysis was presented by Alexander et al. (2003; hereafter A03). 
For this work we used the raw data available from the pub- 
lic archives . Details of the data reduction can be found in 
iLaird et alJ (|200^. The 200 ks Chandra ACIS-I survey of 
the GWS is currently the third deepest blank-field , extra- 
galactic Chandra survey field, after the HDF-N and Chan- 
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Table 1. Fields included in analysis. Col.(l): LBG field name; Col. (2): Type of Chandra observation; Cols. (3, 4): Nominal Right Ascension 
and Declination of Chandra pointings; Col. (5): Galactic column density; Col. (6): Exposure time after GTI, background filtering etc; Col. (7): 
Area of overlap between LBG survey and Chandra observations; Col. (8) Number of LBGs in combined X-ray-LBG survey area; Col. (9) 
Number of spectroscopically identified LBGs with 2>2. 



Field 


Chandra 


RA 


Dec 


Nh 


Filtered Exposure 


Survey area 


No. of 


No. of LBGs 




observation 


{J2000) 


{J2000) 


(lO^^cm-^) 


(ks) 


(arcminutes^) 


LBGs 


with Zspcc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HDF-N 


ACIS-I 


12:36:47.59 


+62:14:08.06 


1.6" 


1862.9 


149.1 


295 


93 


GWS 


ACIS-I 


14:17:43.04 


+52:28:25.20 


1.3* 


190.6 


239.0 


334 


200 


Q1422+2309 


ACIS-S 


14:24:35.61 


+22:55:43.76 


2.7*" 


28.4 


78.1 


292 


103 


SSA22 


ACIS-S 


22:17:28.24 


+00:15:09.59 


4.7*" 


77.8 


105.5 


184 


86 


B20902+34 


ACIS-S 


09:05:33.28 


+34:09:07.83 


2.3'' 


9.78 


41.8 


76 


38 


3C 324 


ACIS-S 


15:49:46.41 


+21:25:19.58 


4.3'' 


38.5 


36.9 


45 


9 



" IStark et alJll992l 

'' iDickev fc Lockmanlll990l 



dra Deep Field-South (CDF-S). A description of the data 
and reduction, inclu dinR a point source list, is detailed in 
iNandra et all Jiooi). 

The remaining LBG fields - Q1422+2309, SSA22, 
B20902+34 and 3C 324 - were each observed once with the 
Chandra ACIS-S instrument at prime focus, with exposure 
times ranging from ~10 ks to ~80 ks. Details of the obser- 
vations are summarised in Tables and |5| The data were 
reduced using the Chandra X-ray Center (CXC) Chandra 
Interactive Analysis of Observations (CIAO) software, ver- 
sion 3.1, and the Chandra calibration database (CALDB) 
version 2.27. The reduction procedure was very similar to 
that employed with the HDF-N and GWS data, but more 
straight forward due to each target being observed just once. 
Images and event files were created for the same four analysis 
bands as used in the HDF-N and GWS analyses: 0.5-2 keV 
(soft band), 0.5-7 keV (full), 2-7 keV (hard) and 4-7 keV 
(ultra-hard). Exposure maps for the LBG survey areas of the 
fields were made using MERGE_ALL with the same represen- 
tative energies that were used for the HDF-N and GWS data 
- namely 1 keV, 2.5 keV, 4 keV and 5.5 keV for the soft, full, 
hard and ultra-hard bands respectively. Filtering to exclude 
periods of unstable background using ANALYZE_LTCRV re- 
moved 0.4 ks, 1.1 ks, 0.1 ks and 4.2 ks from Q1422+2309, 
SSA22, B20902+34 and 3C 324, respectively. 

The final filtered exposure times and aim points of the 
Chandra data are shown in Table Figure shows the 
cumulative survey solid angle as a function of effective ex- 
posure over the entire Chandra /\jQG survey area. 



2.3 Source detection, photometry and matching 

Source detection and photometry was carried out using our 
own procedure which is described in Nandra ct al. (2005). 
As this work is only interested in the X-ray counterparts to 
the LBGs we set the detection threshold to the low value 
of 10"*. This is a less stringent value than one which would 
be acceptable for analysis of the general X-ray source popu- 
lation of our fields. Source detection was carried out in the 
four bands described in §2.2 and a band merged catalogue 
produced. A detection with Poisson probability less than 
10~* is required in at least one band. For all the sources 
in the band merged catalogue accurate source photometry 



Table 2. Table of ACIS-S observations of LBG fields. Each field 
was observed only once. Col.(l): LBG field name; Col. (2): Chan- 
dra observation identification number; Col. (3): Date of beginning 
of observation; Col. (4): Nominal roll angle of satellite (degrees 
East of North). 



Field 




Obs. ID 


Date 


Roll 


Name 






(UT) 


Angle 


(1) 




(2) 


(3) 


(4) 


Q1422+ 


■2309 


367 


2000-06-01 


218.6 


SSA22 




1694 


2001-07-10 


120.3 


B20902- 


+34 


1596 


2000-10-26 


71.0 


30 324 




326 


2000-06-25 


221.1 



is determined using 90 per cent EER apertures in the full, 
hard and ultra-hard bands and 95 per cent EER apertures 
in the soft band. Upper and lower confidence regions for the 
background subtracted counts are calculated for la statist i- 
cal errors (according to equations 7 and 14 of lGehre'l3ll986ll . 
If the Poisson probability in a band is greater than or equal 
to 1.3 X 10""^ (equivalent to a Gaussian Scr detection) then 
we calculate the upper limit to the counts based on that 
same probability. 

The PSFs were calculated at the representative ener- 
gies of each of t he four bands , acco rding to the proce- 
dure described in iNandra et al.l (|200^. In the HDF-N we 
use ex posure weighted averages of the individual PSFs (see 
iLaird ct al. 200 J. To convert count rates to fiuxes we as- 
sume a power- law source spectrum with photon index r=1.4 
and Galactic column density (Table 0. Fluxes in the full, 
hard and ultra-hard bands were extrapolated to the stan- 
dard upper limit of 10 keV. The effects of the Chandra ACIS 
quantum efficiency (QE) degradation (Marshall ot al. 2004)) 
have been taken into account in all fiuxes quoted in this pa- 
per. 

The X-ray catalogues were cross-correlated to the LBG 
candidates using the following procedure. For each field we 
first matched the catalogues using a search radius of 2.0 arc- 
sec to identify possible counterparts. Any astrometric offsets 
between the X-ray and LBG reference frames were then iden- 
tified and the LBG positions corrected accordingly. Astro- 
metric offsets were present in all the fields, with the neces- 
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sary shifts varying from 0.16 arcsec (HDF-N) to 2.06 arcsec 
(GWS). After correcting for the overall offsets the Chandra 
and LBGs catalogues were re-matched using a radius of 1.5 
arcsec. 

The false detection and match rate was assessed by cre- 
ating false LBG catalogues for each field. Using 200 000 
random 'galaxy' positions and matching to the X-ray cata- 
logues using a radius of 1.5 arcsec we find that the probabil- 
ity of detecting an X-ray source at the position of a known 
LBG at random is 0.004, 0.001, 0.002 and 0.001 for the 
HDF-N, GWS, SSA22 and Q1422+2309 fields, respectively. 
Given the number of LBGs and X-ray sources in the fields 
this corresponds to an expected 0.92, 0.42, 0.34 and 0.33 
false matches in the HDF-N, GWS, SSA22 and Q 1422-^2309 
fields, respectively. 



3 BASIC PROPERTIES OF X-RAY 
DETECTED LBGS 

A total of 24 LBGs in the six fields were found to have 
significant X-ray emission - ten in the HDF-N, six in the 
GWS, three in Q 1422-^2309 and five in SSA22a. We did not 
find any X-ray counterparts to the LBGs in fields SSA22b, 
B20902+34 and 3C 324. Twenty of the LBGs have confirmed 
redshifts with 2.211 < z < 3.630. All of the LBGs that were 
identified as QSO or AGN by their rest-frame UV spectra 
in the HDF-N and SSA22 fields have been detected in the 
Chandra data. Table IHTl presents the basic X-ray and optical 
properties of the sources, including coordinates, magni- 
tude, redshift. X-ray fluxes and hardness ratio (HR). 

In the HDF-N, only one of the five least luminous X-ray 
detected LBGs, CXO J123622. 5-^621306 (=HDF-C14), has 
already been reported as an X-ray source in the 2 Ms CDF- 
N catalogue of A03 (source 133). An IR-optic al counterpart 
was d etected in bands HK' through U by iBareer et alJ 
i2003h and a spectrum was taken but no identification 
or spe ctroscopic redshift could be obtained. ISteidel et alJ 
l|2003r classified this LBG as a normal galaxy at z=2.981. 
The remaining four LBGs, CXO J123618.4-F621139 
(=HDF-D7), CXO J1236704.2-H624446 (=HDF-oMD49), 
CXO J123645.0-f621653 (=HDF-M35) and CXO 
J123651. 5+621041 (=HDF-M9), were not in the A03 
catalogue which covered the whole of the CDF-N area and 
had a more conservative detection threshold. As discussed 
in §2.3, less stringent detection thresholds, which result in 
lower flux limits, can be used when searching only for the 
counterparts to limited numbers of known sources as is the 
case in this work. 

The 0.5-2 keV flux of the LBGs varies over two or- 
ders of magnitude from < 3 x 10~^^ erg cm~^ s~^ to 
7.4 X 10~^^ erg cm~^ s"^. The band flux (correspond- 
ing to rest frame ~1800 A emission) covers as large a range, 
from the bright i^=20.48 QSO, HDF-MD39, to very faint 
^=25.27 galaxy, SSA22a-M14, associated with the Lya neb- 
ula. 

Figure |21 shows R magnitude versus 0.5-2 keV flux 
for the detected LBGs. ^% magnitude s were converted to 
Kron- Cousins R using the conversion in lSteidel fc HamiltonI 
l)l993h . Lines of constant X-ray to optical flux are also 
shown , according to the relation from iHornschemeier et aD 
l|200ltl : 
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Figure 2. Soft band X-ray flux vs. R magnitude for the 24 X- 
ray detected LBGs. Black squares and circles denote optically 
identified broad line QSOs and optically identified narrow line 
AGN, respectively. Black triangles denote those LBGs with opti- 
cal galaxy classifications (GAL). Black inverted triangles denotes 
LBGs without optical spectra. The red stars identify HDF-M9 
and HDF-M35 which are GAL LBGs with X-ray fiuxes and spec- 
tra that do not clearly indicate the presence of an AGN and are 
consistent with emission from star formation. The red squares 
show the soft band stacking results for the undetected LBGs in 
the CDF-N and GWS fields. Arrows signify the 3(t upper limit of 
sources not detected in the soft band. ^ band magnitudes have 
been converted to Kron-Cousins R. The diagonal lines indicate 
constant X-ray to optical flux ratios, as indicated. 



log(^^) = log /x + 5.50+^. (1) 

The ratio of X-ray to optical flux is a commonly used tech- 
nique for determining the nature of X-ray sources (e.g. 
iMa ccacaro et alJ Il988l: IStocke et aP Il99ll : ISchmidt etall 
199^" 

Over several decades in flux classic AGN, in- 
cluding luminous narrow line AGN (NLAGN) as well 
as broad line AGN (BLAGN), exhibit X-ray to optical 
flux ratios of — 1 < log (fx/fopt) < -|-1 in both soft and 
hard bands (e.g. ISchmidt et al.lll99i: lAkivama et al.ll2000l: 
iLehmann et alJ l200lt [Alexander et al.ll2001^ ~ At lower X- 
ray to optical flux ratios, mainly in the soft band at 
/o.5-2keV % 10^^ erg s~^ cm~'^, starburst galaxies, nor- 
mal galaxies, low luminosity and heavily absorbed AGN 
emerge (e.g. jGiacconi et all I2OOII: lAlexander et alJ I2OOII : 
iBarger et all 12002: Ba rger et alJ 120031) . Starburst galaxies 
and AGN populate the —2 < log (fx/fopt ) < — 1 range 
(e.g. iBauer et al.ll2002l : lAlexander et alJl200^ while normal 
galaxies which hav e very weak X-ray emission dominate at 
log (fx/fopt) < -2 iHornschemeier et aljliool . 

The X-ray to optical flux of the majority of the LBGs 
falls within the region populated by classic, luminous AGN 
(shaded region). Two of the QSOs and one AGN are sub- 
luminous in X-ray compared to their optical emission for 
typical active galaxies. Half of the LBGs identified as nor- 
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mal galaxies based on their UV spectra have X-ray to opti- 
cal flux ratios typical of NLAGN and BLAGN. Two LBGs, 
HDF-M9 and HDF-M35, exhibit ratios typical of starburst 
galaxies (see §6 for further discussion). As a cautionary note, 
the location of the LBGs in Figure |5| could be affected by 
K-corrections which have been shown to have a redshift de- 
pendent effect on lo g (fx/fopt) values for Seyfert galaxies 
ilPeterson et~ai]l2006l) . 



4 X-RAY SPECTRAL ANALYSIS 

4.1 Hardness Ratios 

To assess the X-ray spectral slope of the LBGs, based on 
broad band X-ray photometry, the effective photon index 
(r) was calculated for each LBG using the observed soft-to- 
hard band ratio. The calculations were made using PIMMS 
v3.6 and assuming only Galactic absorption. The results 
for sources with meaningful soft-to-hard band ratios (i.e. 
excluding sources with only full band detections or those 
with no hard band detections and soft band fluxes close to 
the detection limit) are shown in Table |21 As an alterna- 
tive measure of the X-ray spectrum, the intrinsic column 
density required to produce the observed band ratio for an 
underlying assumed power-law spectrum with F = 2 was 
also calculated (Table 1^1. A F ~ 2 powe r-law is the canon- 
ical spectrum of unobscured AGN (e.g. iNandra fc Pound3 
Il994 : George et al. jOOI3} and additionally well represents 
the average 2-10 keV spectrum of star forming galaxies 
JPtak et al.lll999ll . 

The effective F of the LBGs covers a wide range from 
>2.2 for SSA22a-D13 (indicating a soft, unabsorbed spec- 
trum) to <0.8 for SSA22a-M14 (representative of a flat or 
absorbed spectrum). Alternatively, assuming an underlying 
continuum with a F=2 power-law, the LBGs have intrinsic 
column densities of Nh ~ 0.6 - 3.0 x 10^^ cm"^. 

4.2 Spectral fitting 

To investigate further the effects of absorption, and see if 
there is a correlation between X-ray spectral properties and 
the optical spectral classifications. X-ray spectral analysis 
was carried out. CIAO v3.2.1 and CALDB vS.O.l were used 
for the spectral extraction and XSPEC vll.3.1 software used 
for the spectral analysis. The energy range of all of the anal- 
ysis was restricted to 0.5-7 keV. 

For sources in the GWS, SSA22 and Q1422 fields, the 
source and background spectra were extracted using the 
CIAO tool PSEXTRACT. The source spectra were extracted 
from a circular aperture with a 95 per cent EER calculated 
at 2.5 keV. Local background regions were manually selected 
to avoid contamination by nearby X-ray sources. Typically, 
the background was extracted from three to four large re- 
gions surrounding the source. 

As a result of the different aim points, roll angles and 
observing modes of the HDF-N observations the spectra in 
this field could not be extracted using PSEXTRACT in one 
step. Instead, source and local background regions were de- 
fined as before and spectra were extracted for each of the in- 
dividual observations. The spectra were then summed using 
the standard FTOOLS jBlackburnll995fl routine MATHPHA. 



Table 3. Effective photon indices (F) and, alternatively, effec- 
tive intrinsic column densities (Njj) calculated from the hard-to- 
soft band ratios, for sources with meaningful limits. Col.(l): LBG 
name; Col. (2): Effective F over 0.5-10 keV assuming only Galac- 
tic absorption. Col. (3): Intrinsic column density calculated for an 
underlying power-law spectrum with F = 2 in units of 10^"^ cm"-^. 



LBG 


F 


Nh (for F = 


= 2) 






(10^3 cm- 






{2") 


('3') 




HDF-D7 


< 1.0 


> 2.7 




HDF-C14 


1 n+0.3 

^•'-'-0.2 


2 O+l-l 




HDF-OC34 




1 q+0.5 
1-^-0.5 




HDF-CIO 




1-^-0.3 




HDF-MD34 


1 7+0.1 


-0.3 




HDF-MD12 


1 2+°-l 


r, O + 0.3 

^-'-'-0.2 




HDF-MD39 




6+°-^ 
^-°-o.i 




Q1422-C73 


1 O+0.4 
^•■^-0.3 


9 c:+2.5 




SSA22a-D13 


> 2.2 


No abs. 




SSA22a-D12 


1 3+0.4 
^•■^-0.3 


9 q+1.9 
^-■^-1.5 




SSA22a-M14 


< 0.8 


> 5.6 




GWS-MD106 


l-°-0.1 


o.8tr4 




GWS-D54 


1 4+0-2 






GWS-M47 


1 1+0.3 
1-1-0.2 






GWS-C50 


1-^-0.3 







The instrument response and matrix files were coadded us- 
ing the FTOOLS routines ADDARF and ADDRMF, weighted 
by exposure time. The response files for the first three HDF- 
N observations (IDs 580, 966 and 967) are different to those 
for the remaining observation IDs as a result of the higher 
focal plane temperature for these observations. Therefore 
when combining the individual spectra only the 17 obser- 
vations taken at — 120°C were included, reducing the total 
exposure of the spectra by 161.7 ks. 

Most of the sources have limited counting statistics 
(< 200 counts) so the spectr al analysis was performed using 
the C-statistic (|Cas3[l97^, which was specifically devel- 
oped to extract information from spectra with low numbers 
of counts. The sample was limited to sources that had a 
minimum of 10 total counts in the extracted spectra, result- 
ing in the exclusion of SSA22a-M8, SSA22a-M14, SSA22a- 
MD14 and Q1422+2309b from the spectral analysis. AU 
of the spectra analysed using the C-statistic were grouped 
into fixed width bins of 4 channels (~ 700 eV), which im- 
proves the processing times in XSPEC. For those sources 
with > 100 counts in their spectra (HDF-oC34, HDF-MD34, 
HDF-MD12, HDF-MD39 and GWS-MD106) standard 
spectral fitting was also performed. In this case the spec- 
tra were grouped to have a minimum of twenty counts per 
bin, required for approximately Gaussian statistics. 

The LBG spectra were also searched for the presence of 
iron Kq line emission by looking for evidence of an excess of 
counts at rest frame 6.4 keV, as compared to that expected 
from the best fitting model for each source. None of the 
LBGs showed significant evidence of Fe K line emission. 
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Table 4. X-ray spectral fits using the C-statistic iCashlll97fli^ for sources with greater than 10 counts. All errors correspond to the 
90 per cent confidence level. Col.(l): Optical spectral classification; Col. (2): LBG name; Col. (3): Number of counts in spectral fitting; 
Col. (4): Eff'ectivc F for fits assuming only Galactic absorption; Col. (5): Intrinsic column density, Nh, in units of 10^'' cm~^, for fits 
assuming a fixed intrinsic F = 2 power-law spectrum with Galactic absorption; Col. (6): Improvement in the C-statistic for an intrinsic 
F = 2 spectrum by allowing the absorption shown in Col. 4 compared to zero intrinsic absorption; Col. (7): Intrinsic column density 
found from simultaneous fitting of whole sample, where the best fitting photon index was F = 1.96. For HDF-MD39 and HDF-MD12 
(which were excluded from the simultaneous fitting) the results are for a fixed F = 1.96 power-law. 
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4.3 C-statistic and fitting of individual LBGs 

For each of the LBGs the data were fitted in several ways 
using the C-statistic to assess the level of intrinsic absorp- 
tion and the effective power-law photon index, F. First, to 
determine the effective photon index F of the spectra, the 
data were fitted with a power-law model with Galactic ab- 
sorption using the Nh values in Table Q and zero intrinsic 
absorption (column 4 of Table 2J- Secondly, the data were 
fitted with a power-law model absorbed by both an intrinsic 
column density at the source redshift and a Galactic col- 
umn density (Table |11 column 5). In this case the photon 
index was fixed to F=2 in order to constrain the intrinsic 
Nh required to produce the observed spectra. Adopting a 
fixed r=2 spectrum (the canonical spectrum of unobscured 
AGN) is an approximation used to assess the intrinsic col- 
umn density. For sources with a fiatter intrinsic spectrum 
this will result in an overestimation of Nh. The errors in 
the fitted values of both F and Nh correspond to the 90 per 
cent confidence level for 1 interesting parameter (AC= 2.7). 
As a further measure of the intrinsic absorption, the data 
were also fitted with fixed r=2 with only Galactic absorp- 
tion and the improvement in the C-statistic between this fit 
and the previous fit allowing intrinsic absorption was calcu- 
lated (Table |1] column 6). Large values of AC-stat indicate 
that the inclusion of an intrinsic column density produced 
substantial improvement to the fits. 

In addition to the individual C-statistic fits, simultane- 



ous fitting of the sample was also performed to determine 
the mean photon index of the LBGs as a whole. In this fit 
both the photon index and intrinsic column density were 
free parameters, but F was fixed to be the same for all the 
LBGs. HDF-MD39 and HDF-MD12 were excluded from the 
simultaneous fitting to prevent the signal being dominated 
by these two very bright LBGs. The best fitting photon in- 
dex for the LBG sample was found to be F = IM^o^^l- 
The intrinsic column density of each LBG obtained from 
the spectral fitting is given in Table^] column 7. The values 
of Nh found from the simultaneous fitting are very simi- 
lar to those that were found when assuming a fixed F = 2 
spectrum. 

The spectra of the five LBGs with >100 counts were 
fitted using the statistic. Each were fitted with a fixed 
F = 2 power-law spectrum with both intrinsic and Galactic 
absorption, as was carried out with the C-statistic. The de- 
rived Nh values along with the reduced s-id probability 
are shown in Table |K1 columns 2-4. The three LBGs with 
>200 counts were also fitted with F as a free parameter. 
The results are shown in Table |S] columns 5-8. Strong con- 
straints were able to be placed on the values of F and Nh 
for HDF-MD39 (Figure 1^1 and HDF-MD12 (Figure Hi. 

The results of the various spectral fits show that the 
majority (~ 60 per cent) of the LBGs have spectra con- 
sistent with an unobscured AGN, although several sources 
(particularly those classed as galaxies optically) are too faint 
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Table 5. X-ray spectral fits using the x-^-statistic. For all sources, fits were performed assuming Galactic absorption and a fixed 
intrinsic F = 2.0 power-law spectrum. For HDF-oC34, HDF-MD39 and HDF-MD12, additional fits were performed where F was 
also a free parameter. Col.(l): LBG name; Col. (2): Fitted intrinsic column density, Njj, in units of 10^^ cm~^. Uncertainties refer to 
Ax^ = 2.71, corresponding to the 90 per cent confidence level for one interesting parameter; Col. (3): Reduced x^; Col. (4) x^ probability; 
Col. (5): Fitted photon index, F; Col. (6): Fitted intrinsic column density (units of 10'^^ cm"'^). The uncertainties in Cols. 5 and 6 refer 
to Ax^ = 4.61, corresponding to the 90 per cent confidence level for two interesting parameters; Col. (7): Reduced x^; Col. (8) 
probability. 
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Figure 3. The X-ray spectrum of CXO J123622.9-I-621526 
(=HDF-MD39), binned to 20 counts per bin. The fitted model 
of Galactic absorbed power-law emission with intrinsic F = 
l.SO^QQy and negligible intrinsic absorption is shown. Data-to- 
model residuals are shown in the bottom panel in units of cr. The 
spectrum is consistent with no intrinsic absorption. 



Figure 4. The X-ray spectrum of CXO J123719.8-I-620955 
(=HDF-MD12), with a minimum of 20 counts per bin. The fit- 
ted model of Galactic absorbed power-law emission with intrinsic 
F = 1.92t[5;^4 and intrinsic column density Nh = l-65to gg ^ 10^^ 
is shown (90 per cent confidence range for two interesting param- 
eters). The bottom panel shows the the data-to- model residuals 
in units of cr. 



to be able to place reliable constraints. Seven LBGs show 
evidence of significant obscuration (SSA22a-D13, SSA22a- 
D12, HDF-MD12, Q1422-MD109, Q1422-C73, GWS-M47 
and HDF-D7). Of the seven LBGs classified as broad line 
QSOs, two show evidence of intrinsic absorption (SSA22a- 
D12 and -D13). The brightest LBG in the sample, the QSO 
HDF-MD39, is well fit by a F = 1.80 power-law with no 
absorption (Table |^ and Figure Ol. Four of the five LBGs 
classified optically as narrow line AGN exhibit evidence for 
significant intrinsic absorption, including HDF-MD12 (Fig- 
ure HJ. The remaining AGN classified LBG, HDF-oMD49, 
may also be obscured but is too faint to provide reliable 
constraints. Longer wavelength data lend support to and 
obscured hypothesis - the UV spectrum shows evidence of 
self-absorption, the source has a very bright 24 micron flux 
(~300 ^Jy) and a power-law SED in the IR (C. C. Steidel, 
priv. communication) . 



4.4 Simultaneous fitting of LBG subsets 

To examine the average spectral properties of the LBGs 
as a function of optical spectral classification, simultane- 
ous fitting was carried out for separately for objects classi- 
fied as broad or narrow line AGN, and galaxies. For each 
class two simultaneous fits were performed: the first with a 
fixed F = 1.96 power-law spectrum with intrinsic absorp- 
tion, and the second with a power-law spectrum with in- 
trinsic absorption, allowing F to be a free parameter. Again, 
the two brightest LBGs were excluded from the fits to pre- 
vent the signal being dominated by these individual sources. 
The results are shown in Table 1^1 The mean spectrum of 
the QSOs is consistent with a F ~ 1.9 power-law with at 
most a small amount of absorption (1.9l^ g x 10^^ cm~^). 
The mean spectrum of the AGN shows more absorption 
(Nh = l-Tli'a ^ 10^^ for a F = 1.8 power-law), consistent 
with the results of the individual LBGs in §4.3. The mean 
spectrum of the galaxy class is poorly constrained. 
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Table 6. Simultaneous C-statistic fitting for LBG optical classes. 
HDF-MD39, HDF-MD12 and sources with <10 counts in the 
spectrum are excluded. Col.(l): Optical spectral classification; 
Col. (2): Intrinsic column density for fits with a fixed F = 1.96 
power-law spectrum; Cols. (3 and 4): photon index and intrinsic 
column density Njj (units of 10^"^ cm~-^) for two free parameter 
fits. All errors correspond to the 90 per cent confidence level. 
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5 STACKING OF UNDETECTED LEGS 

The mean X-ray properties of the LBGs too weak to be 
directly detected are determined by employing a stacking 
technique, shown to be a useful and successful tool with 
Chan dra data (e.g. Brandt et al. 2001; Hornschcmoicr ct al. 
l200ll : iLaird et alTl2005f) . Previous stacking analysis of 2 ~ 3 
LBGs in the HDF-N, using the 1 Ms data (N 02) and a dif- 
feren t sample of galaxies in the 2 Ms data jLehmer et alJ 
I2OO5I) . have yielded highly significant detections. The emis- 
sion from these X-ray weaker LBGs is thought to be dom- 
inated by star formation process instead of AGN and the 
results have, for instance, being used to calculate X-ray de- 
rived SFRs. The X-ray derived SFRs from stacking have 
been shown to be consiste nt with rad io and extinction cor- 
rected UV estimates (e.g. iReddv fc Steidel 2004') . Here we 
perform stacking analyses of the undetected LBGs in the 
HDF-N and the GWS. 

5.1 Stacking procedure 

The stackin g procedure used in this work is identical to that 
described in lLaird et alJ l|20 05') and we include only a brief 
outline here. In both the HDF-N and GWS fields all LBGs 
with an X-ray counterpart are excluded from the stacking, as 
are those with an unassociated close by X-ray source. Source 
counts are extracted from the optical LBG positions (cor- 
rected to account for the astrometric offset between the X- 
ray and optical reference frames, §2.3) using a circular aper- 
ture and are summed to find the total counts for the LBGs in 
the sample. The background counts are estimated using two 
methods. First, the LBG positions are randomly shuffled by 
5-10 arcsec and background counts are extracted from an X- 
ray source-masked image. Secondly, background counts are 
extracted from random positions anywhere over the field of 
view (excluding areas with no exposure or markedly differ- 
ent exposure from that of the LBG positions, in the case of 
the HDF-N). This is repeated 1000 times for each LBG posi- 
tion. The background counts are then summed and scaled to 
the same area as the source extraction to find the net source 
counts. To identify the X-ray sources for the source-masked 
images we repeated the source detection procedure described 
above (§2.3) using a more stringent probability threshold of 
10~^, for which the numbers of spurious sources over the 
entirety of both survey areas are expected to be small. 

The size of the extraction aperture used, as well as 



the radius from the Chandra aim point within which to in- 
clude LBG positions, affects the strength and accuracy of the 
stacking signal. We adopt an empirical approach to deter- 
mining the optimal extraction and inclusion radii by testing 
10 fixed extraction radii between 0.75 and 3.0 arcsec and 

6 inclusion radii between 5 and 10 arcmin off- axis and se- 
lecting the radii yielding the maximum signal-to-noise ratio. 
In the HDF-N we chose an extraction radius of 1.5 arcsec^ 
and an inclusion radius of 9 arcmin (Figure |^. In the GWS 
an extraction radius of 1.5 arcsec and inclusion radius of 

7 arcmin produced the strongest signal-to-noise ratio (Fig- 
ure 0. In both fields the shuffled and random background 
methods produced similar results and in this paper we quote 
all results using the shuffled positions, which should better 
account for local variations in the background level. 

In converting the stacking count-rates to fluxes a power- 
law spectrum with F = 2.0 and Galactic absorption was 
assumed, ap propriate for sou rces with high star formation 
activity (e.g. lPtak et al.lll99<ih . 



5.2 Stacking Results 

The results of the stacking in the HDF-N and GWS are 
shown in Table Q Stacking the soft band emission from 
the 277 undetected LBGs within 9 arcmin from the Chan- 
dra aimpoint produces a highly significant detection with 
a signal-to-noise ratio of 5.9 (where signal-to-noise ratio is 
defined as S/\/S -|- B and S and B are the net source counts 
and background counts, respectively). An average of 0.58 net 
counts per galaxy are detected, with 2.51 mean background 
counts in each extraction cell. Figure |7[a) shows the distri- 
bution of total counts (S-fB) in each LBG extraction cell. 
The distribution ranges from zero to ten counts per LBG 
and is not dominated by a few sources, but rather is well 
representative of the full sample. The mean flux per HDF-N 
LBG is 2.36±0.43x 10"^^ erg s'^cm'^ (0.5-2 keV). Stacking 
the soft band emission from 226 LBGs in the GWS also pro- 
duced a detection, albeit a marginal one. An average of 0.12 
net counts per LBG were detected, with 0.31 background 
counts in each detection cell, resulting in a signal-to-noise 
ratio of 2.81. The total count in cell distribution for the GWS 
is shown in Figure|7[b) and covers a much smaller range than 
the HDF-N, as would be expected given the shorter expo- 
sure time. The mean soft band flux per LBG was found to 
be 4.34 ± 1.54 x 10"^* erg s'^cm"^ (0.5-2 keV). The mean 
flux in the GWS is 84 per cent larger than in the HDF-N, 
but consistent within the errors at the 1.2a level. However, 
a raw comparison of the two results is misleading because 
of the different flux limits in the GWS and HDF-N. In or- 
der to provide a fairer comparison between the two samples 
the HDF-N stacking was repeated with the inclusion of the 
six detected LBGs with fluxes below the flux limit in the 
GWS. The mean flux per LBG in this sample was found to 
be 3.23 ±0.44 x 10"^** erg s"^cm"2 (0.5-2 keV), completely 
consistent with the GWS result. 



^ We note the N02 contains an error regarding the optimal ex- 
traction radius used. The paper states that the optimal extraction 
radius found was 2.5 arcsec; in fact, 2.5 arcsec was the optimal 
extraction diameter. 
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Figure 5. The signal-to-noise ratio vs. (a) extraction radius in arcseconds and (b) inclusion radius in arcminutes for the LBGs in the 
HDF-N. The results are shown for two different background methods - shuffled (circles) and random (squares) positions. The signal- 
to- noise ratio is an inverse measure of the fractional error on the flux and is given by (S/\/S + B), where S and B are the net source 
counts and background counts respectively. The vertical dashed lines denote our chosen extraction and inclusion radii of 1.5 arcsec and 
9 arcmin, respectively. 
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Figure 6. The signal-to-noise ratio vs. (a) extraction radius in arcseconds and (b) inclusion radius in arcminutes for the LBGs in the 
GWS. The symbols and axes are the same as for Figure |S] The vertical dashed lines denote the chosen extraction and inclusion radii of 
1.5 arcsec and 7 arc, respectively. 



Stacking the hard band emission did not result in a 
significant detection in either field. 

5.3 X-ray and UV correlations for stacked LBGs 

The soft band stacking signal in the HDF-N is sufficiently 
strong to allow the sample to be split into bins according to 
optical properties and the mean X-ray properties determined 
for each bin. This has been shown to be very effective in the 
HDF-N for BBGs at z ~ 1, allowing correlations between X- 
ray a nd rest-frame UV emission to be examined iLaird et alJ 
12005^ . In order to achieve reasonable signal-to-noise ratios 
the LBG sample was split into three bins of approximately 
90 galaxies according to magnitude, which corresponds 
to rest-frame ~ 1800 A emission at z ~ 3. The stacking was 
performed according to the procedure described above, using 
a 1.5 arcsec extraction radius and only including galaxies 
within 9 arcmin. 

The subset stacking results are shown in Table Q a^nd 
Figure |H| A significant detection was found for the brightest 
^ bin and marginally significant detections found for the 
remaining bins. In each case the fiux was sufficiently con- 
strained to allow an assessment of any correlations. As can 
be seen from Figure |H| within the limited statistics and dy- 
namic range of this LBG sample, we find no evidence for a 
correlation between ^ magnitude and soft X-ray flux, cor- 
responding to rest-frame 1800 A and 2-8 keV emission. 



6 DISCUSSION 

We have presented results from Chandra X-ray observa- 
tions of 6 fields which have been surveyed with deep optical 
photometric and spectroscopic observa tions designed to se - 
lect Lyman break galaxies at z ~ 3 dSteidel et al.ll20o3) . 
These observations have approximately doubled the num- 
ber of known X-ray d etect ions of LBGs to a total of 24 
(c.f.lBrandt et a l.'200l': N02: lNandra. Laird, fc Steidei2005l : 
iLehmer et al J 12005). The raw fraction of X-ray detected 
LBGs in the entire sample is 2 per cent, but it should be 
borne in mind that many of the Chandra observations pre- 
sented here are very shallow, and not sufficiently sensitive to 
probe deep into the X-ray luminosity function. Taking the 
HDF-N alone, which has the deepest data, we estimate an 
X-ray detection fraction in LBGs of 3 ± 1 per cent, similar 
to that found bv lSteidel et alJ J2002h . 

It is therefore immediately interesting to note that the 
fraction of LBGs hosting a n AGN is much lower t han that of 
submm selected galaxies l|Alexander et alJl2005h . in which 
more than half the sources are likely to harbour an ac- 
tively accreting black hole. The importance of this result 
is at least twofold. Firstly, there remains considerable un- 
certainty about the contribution of AGN to the bolometric 
luminosit y of FIR select ed galaxies at the highest luminosi- 
ties (e.g. IVeilleux et alj fl995: Farrah ct al. 200^. No such 
uncertainty would appear to be pertinent for the LBGs. Sec- 



10 E. S. Laird et al. 



Table 7. Stacking results of undetected LBGs in the HDF-N and GWS, including results for sub-samples based on rest-frame 
1800 A emission in HDF-N. Col.(l): Field. Col. (2): Galaxy sample. Col. (3): Observed-frame energy band. Col. (4): Number of galaxies 
included in stacking sample, taking into account rejected galaxies as described in §2.5. Col. (5): Mean ^ magnitude. Col. (6): Mean 
redshift. Col. (7): Signal-to-noise ratio [S/\/S~+B], where S and B are the net source and background counts respectively. Col. (8): X-ray 
flux per galaxy in units of 10~^* erg cm~^ s~^; 0.5—2 keV for soft band, 2-10 keV for hard band. Col. (9): X-ray luminosity per galaxy 
in the 2-10 keV band, derived from soft band flux assuming r=2.0 and Galacti c Njj. 10-50 keV luminosity is given for hard band, 
derived using r=1.4. Col.(lO): SFR from 2-10 keV luminosity iRanalli et al.l2003l) . Errors are statistical only. Col. (11): Ratio of X-ray 
derived SFR to UV SFR, uncorrected for attenuation. 
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24.81 


3.00 


7.3 


3.23 ± 0.44 


2.78 ± 0.38 


55.6 ± 7.6 


5.4 ± 0.7 


GWS 


All undetected 


Soft 


226 


24.70 


2.94 


2.81 


4.34 ± 1.54 


3.60 ± 1.28 


72.0 ±25.6 


6.4 ± 2.3 


GWS 


All undetected 


Hard 


229 


24.69 


2.94 


<0.0 











" Includes directly detected LBGs with fluxes below the flux limit of the GWS. 
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Figure 7. Soft band counts distribution for the undetected LBGs 
in (a) the HDF-N and (b) GWS. The vertical dashed line denotes 
the mean background counts per extraction cell, derived via the 
shuffle background method. 



ondly, the results for the submm galaxies have been inter- 
preted as a strong linkage between intense star formation 
and AGN activity at high- redshift iAlexander et al]|2005t 
see also Page et al. 2001. 2004). We find no evidence that 
LBGs, which are certainly galaxies in which active star for- 
mation is occurring, are also preferentially active in nuclear 
black hole accretion. 

The X-ray-to-optical flux ratios of the LBGs cover al- 
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23.5 



24.0 



24.5 



25.0 



25.5 



26.0 



Figure 8. X-ray stacking results for undetected LBGs in the 
HDF-N split into subsets based on ^ band magnitude. At z ~ 3 
.5^ corresponds to rest-frame ~ 1800 A emission. The x-axis error 
bars are the standard deviation of values in a given bin and the y- 
axis error bars are the Poisson errors from the stacked soft band 
counts. There is no correlation between soft X-ray flux and Si 
magnitude. 



most the full range observed in X-ray selected samples. Ob- 
jects with the highest X-ray-to-optical ratio (e.g. EXOs; 
ipXockomocr ot al. 2004) are not present but this is clearly not 
surprising given this is an optically selected sample. Most of 
the X-ray detected LBGs fall within the range expected for 
AGN. A few objects show lower X-ray-to-optical flux ratios, 
but four of these are clearly AGN also, based either on their 
optical spectroscopic properties. X-ray luminosity, or X-ray 
spectral properties. LBGs hosting AGN therefore display a 
range in X-ray to-optical flux ratio of at least two orders of 
magnitude. 

There are two X-ray detected LBGs, HDF-M9 and 
HDF-35 which exhibit very low X-ray-to-optical flux ratios, 
typical of starburst galaxies. Furthermore, both objects have 
soft X-ray spectra, low X-ray luminosity (~ 1 — 3 x 10*^ erg 
s"'^) and have no indication of AGN activity in their op- 
tical spectra. It is possible that their X-ray emission is 
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simply due to low luminosity AGN activity. It also seems 
quite plausible, however, that the X-rays from these two 
objects are powered by intense star forming activity, with 
the emi ssion being predominantly from X-ray binaries and 
hot gas jPavid. Jones, fc Formanll992h . Using the standard 
conversion s between X-ray luminosity and star-formation 
rate (N02 ; iRanaUi et all2003l: iGrimm. Gilfanov. fc SunvaevI 
l2003l) . we infer star formation rates for the two objects of 
292t^;^i M© yr"^ (HDF-M9) and 52lti42 Mq yr"^ (HDF- 
M35). Clearly these objects would then be considered quite 
extreme, wit h SFRs comp arable to bright submm selected 
galaxies (e.g. Ilvis on et al] l 2000') or Hyper luminous IRAS 
galaxies iRowan-Robinsonl l20001 . These high SFRs should 
be compared to the typical, extinction-corrected average val- 
ues for LBGs, which we find in th is work to be ~ 40 — 50 Mq 
yr~^ (see also lSteidel et al.ll999l : N02). Given the large sam- 
ple under consideration, and the fact that we expect a wide 
range of SFRs in the LBGs, it seems quite reasonable that 
there are extreme objects with SFRs a factor of several 
higher than the mean. Subsequent analysis of Spitzer IRAC 
and MIPS 24^im observations yielded SFRs of ~170 M© 
yr"^ (HDF-M9) and ~360 M© yr"^ (HDF-M35) and sup- 
ports a star formation hypothesis for both these LBGs (C. 
C. Steidel, priv. communication). 

Whatever the origin of the X-ray emission in these two 
objects, our analysis has revealed new information about 
the AGN population at high-redshift. We find that the X- 
ray detected LBGs display a wide range of luminosities from 
log Lx = 42 — 45. T his is typ ical of the AGN population at 
lower redshift (e.g. ICowie e t al. 2002; Rosati ct al. 20o3). 
The optical spectral classifications of the X-ray detected ob- 
jects are quite diverse, comprising approximately 1/3 broad 
line AGN, 1/3 narrow line AGN and 1/3 with spectra typical 
of normal star forming LBGs. These proportions are again 
very similar to AG N in X-ray selected samples in general 
jBareer et al.ll2003f) . The implication is that the range of 
properties of AGN in the LBG population at 2; ~ 3 are very 
similar to those at lower redshift. In particular, our results 
support the idea that there is a large population of low lu- 
minosity "Seyfert" level AGN of all optical types at 2 = 3 
JSteidel et alJl2002t iNaiidra. Laird, fc Steidelll2005^ . 

X-ray spectral analysis of the LBGs shows evidence for 
significant absorption in several objects. This may be con- 
sidered surprising given both the limited photon statistics in 
the spectra, and the fact that at 2 ~ 3 we sample the spectra 
only ai E > 2 keV in the rest frame, so that only the heaviest 
absorption is detectable. Formally, we find evidence for ab- 
sorption in about half the sample, but given the above con- 
siderations the true fraction could be higher. Splitting the 
objects according to optical spectroscopic class and fitting 
simultaneously, we find the clearest evidence for absorption 
in the narrow line objects, consistent with standard unifi- 
cation schemes where the s ame material obscures bot h the 
X-rays and the optical BLR llAntonucci fc MilleI^ll985^ . The 
columns are typically 10^^ cm"^^ a litt le lower than thos e 
for type II Seyferts in the local universe jAwaki et al.ll99ll) . 
Once the X-ray fluxes are corrected for absorption (TablelSJ, 
a number of the narrow line objects fall into the luminos- 
ity and obscuration range where they would be classified as 
candidate "type II QSOs". For example, four objects have 
Lx > 10** and A^h > 10^ ^ cm"^. These have been used as 
typical dividing lines fe.g. iMainieri et al.ll2002l : iBrusa et alJ 



Table 8. X-ray luminosities of the LBGs. CoL(l): LBG name; 
Col. (2): Optical classification; Col. (3): Rest- frame 2-10 keV lu- 
minosity, in units of lO*"^ erg s~^; Col(4): Absorption-corrected 
rest-frame 2-10 keV luminosity, for sources with evidence for ab- 
sorption in their spectrum. 



LBG 



(1) 



Optical 
type 
(2) 



L2-IO kcV 
uncorrected 

(3) 



L2-IO kcV 
corrected 

(4) 



HDF- 
HDF- 
HDF- 
HDF- 
HDF- 
HDF- 
HDF- 
HDF- 
HDF- 
HDF- 



D7 

C14 

oC34 

CIO 

MD34 

oMD49 

MD12 

MD39 

M35 

M9 



GAL 
GAL 
QSO 



AGN 
AGN 
QSO 
GAL 
GAL 



0.42 
0.64 
5.05 
1.13 
2.24 
0.67 
6.84 
49.79 
0.17 
0.26 



-0.09 
-0.08 
+0.15 
0.12 
+0.46 
0.42 
+0.19 
0.16 
+0.27 
0.25 
+0.22 
0.17 
+0.44 
0.42 
+ 1.05 
1.00 
+0.07 
0.05 
+0.10 
0.07 



11.45 
1.02^ 



+6.89 
^11.04 
-1.48 



4.43 



+86.63 



16.81 



+ 11.24 



0.44 



+5.92 
1.79 



Q1422-|-2309b 
Q1422-MD109 
Q1422-C73 


QSO 
AGN 
AGN 


1 n 7r:+5.82 
lU. (0_3 gg 

Ol +3.56 
^•■^^-2.58 

17.83t«;« 


19.431^71 
47.641^,^° 


SSA22a-D13 

SSA22a-M8 

SSA22a-D12 

SSA22a-MD14 

SSA22a-M14 


QSO 

QSO 
GAL 
GAL 


14 1Q+3.93 
l^-lo_3.i4 

7.5il?:f5 

J^-»J_0.83 
^•'°-1.02 


28.901^*0°' 
16.58l«;;| 



GWS-MD106 


QSO 


16.78t?-J0 


21-051^:^1 


GWS-D54 


QSO 






GWS-M47 


AGN 


4.88t{fr 


12.50l|:g 


GWS-MIO 




9 0-1+0.88 
^■■^1-0.66 




GWS-0MDI3 


QSO 


4 26+^-2«' 




GWS-C50 


GAL 


2 ■iP.+O-S^ 





I2005li although it should be noted that absorbing columns of 
lO^'^ cm~^ should be considered on the very low side for bona 
fide Seyfert 2s ( Awaki e t al. 1991). The nature and origin of 
the absorption in these lightly-obscured objects, and hence 
in the LBGs, remains uncertain and could easily be related 
to galactic-scale dust an d gas rather than a standard torus 
jMaiolino fc Ri ekc* 1995) . It should further be noted that at 
least two objects show QSO luminosity and heavy obscura- 
tion in excess of 10^'' cm~^, but have broad emission lines in 
their optical spectra. These would be classified as candidate 
type II QSOs based solely on their X-ray properties (e.g. 
IMainieri et alJl2002h but are clearly not type lis by defini- 
tion. How these objects fit into the unification schemes is 
currently unclear, but it is possible either that the obscura- 
tion is purely nuclear, affecting the X-rays only, or that the 
torus is relatively warm and ionized, and hence dust-free. 

Taking a naive flux-luminosity conversion we flnd that 
the narrow line AGN sample is less luminous in the X-rays 
than the broad line QSOs (Table This effect seems to 
be primarily due to the more common incidence of heavy 
absorption in the AGN. After correction for obscuration we 
flnd the average luminosity of the broad and narrow line ob- 
jects to be very similar. Objects classified as galaxies from 
optical spectroscopy (i.e. without AGN signatures) are, on 
the other hand, roughly an order of magnitude less lumi- 
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Table 9. Mean 2-10 keV luminosity of LBG optical classes. 
Col.(l): Optical spectral classification. Col. (2): Mean 2—10 keV 
luminosity of spectral group, uncorrected for X-ray absorption. 
Col. (3): Mean 2-10 keV luminosity of spectral group, corrected 
for X-ray absorption. 



Optical 


Uncorrected L2_iofceV 


Corrected L2_iofceV- 


Class 


(10"'^ erg s-l) 


(lO''^ erg s-l) 


(1) 


(2) 


(3) 


QSO 


15.05 ± 14.69 


18.56 ± 15.03 


AGN 


7.71 ± 6.35 


20.16 ± 16.38 


GAL 


1.22 ± 1.09 


3.39 ±4.06 



nous in the X-ray than optically-identified AGN (Table EJ. 
This conclusion holds even after absorption is accounted 
for. The galaxy-classified objects tend not to show evidence 
for absorption (with the clear exception of HDF-D7), but 
small photon statistics prevent definitive conclusions. It is 
therefore possible that large-scale dust obscured both the 
broad and narrow lines in these objects. It seems more likely 
based on our results, however, that the lack of optical AGN 
signatures is due to their low intrinsic luminosity, so that 
the galaxy light dominates over the AGN in the optical 
jMoran. Filippenko. fc Chornockll2002h 

A final note regarding the X-ray spectral properties of 
the directly detected LBGs is that their mean continuum 
spectral index of F = 1.96 is remarkab ly similar to that of 
local Seyferts llNandra fc Founds! Il99ll . We therefore find 
no evidence for evolution of the continuum-generating mech- 
anism with redshift. 

The mean properties of the LBGs not directly detected 
in the deep HDF-N and GWS observations have been de- 
termined by stacking. This technique has been successfully 
applied in the past to constrain LBG X-ray emission which 
is thought to be dom inated by star-form ing pr ocesses, rather 
than AGN activity llBrandt et alJl200ll: N02; iLehmer etal] 
I2OO5I) . The other data sets in this study are too shallow to 
constrain the properties of such sources and, indeed even the 
200 ks GWS observation is of limited value. Analysis in the 
deeper HDF-N field shows a population of LBGs, the ma- 
jority of which are AGN, with X-ray fluxes just below the 
flux limit of the GWS observation. When these sources are 
included in the HDF-N stacking sample they dominate the 
average stacked fiux. This shows that there is a clear risk of 
contamination in the GWS signal by X-ray sources that are 
just below the flux limit for direct detection, and whose emis- 
sion is likely dominated by black hole accretion rather than 
star formation. The best estimate of the stacked flux in the 
GWS, while poorly constrained, is about 80 per cent higher 
than that of the HDF-N, probably due to this contamina- 
tion. This shows the limited usefulness of the 200 ks GWS 
observations in determining the properties of star-forming 
LBGs and henceforth we restrict our discussion to the HDF- 
N. 

We find a highly significant signal in the soft band when 
stacking all 277 undetected LBGs in the HDF-N. There is 
no detection in the hard band. This is primarily due to the 
lower ACIS sensitivity in the hard band, but it shows at least 
that the stacked objects are consistent with a relatively soft 



X-ray spectrum expected from star formation. The inferred 
flux per object is entirely consistent with previous stack- 
ing analysis, and specifically that inferred by N02 using the 

1 Ms data. This is important because it shows that the sig- 
nal detected by N02 was not dominated by a few bright, 
AGN-dominated sources. It should be noted, however, that 
at least one clear X-ray AGN (HDF-C14) has been detected 
in the 2 Ms data that was not detected in 1 Ms. The case 
of HDF-0MD49 is particularly inte resting because this ob- 
ject has already been identified by ISteidel et alJ ll2002ll as 
an AGN based on optical spectroscopic signatures, but was 
not detected in the 1 Ms Chandra data. While we have now 
found a detection in the 2 Ms observation this highlights the 
very wide range of X-ray-to-optical flux ratio in the LBG 
AGN discussed above. 

The mean star formation rate inferred from the X-rays 
for the LBGs is found to be 42 ± 8 Moyr"^ Using the 
ratio of the X-ray derived SFR (which should be roughly 
independent of extinction) to the UV derived SFR (uncor- 
rected for extinction) , the inferred UV extinction correction 
at 1800 A is found to be 4.1 ± 0.8 (Table Q. This is in 
good agreement with the extinction correction inferred pre- 
viously by N02, and with that determined for ga laxies at 

2 ~ 2 jReddv fc SteideJ 12004 iReddv et al.l 120051) . It also 

agrees with, and validates, the extinction corrections em- 

I I J ' I, 

ployed by Stcidel et al. (1999;, which were calculated using 

the Calzetti (1997 ) reddening law with E{B - V) = 0.15, 
typica l for the LBGs. It also agrees well with the analy- 
sis of IPettini et alJ lll99^ . which used H/3 luminosity as 
a comparison SFR measure, and th e UV/X-ray analysis 
ofl Seibert. Heckman. fc MeureJ (|20o3). This extinction cor- 
rection is however on the lower end of an estimate by 
IViih. Witt, fc GordonI (|20o3), who found the UV attenua- 
tion factor for a large sample of 2 ~ 3 LBGs to be luminosity 
dependent and in the range 5.9 to 18.5. 

Our sample is now large enough, and the stacking sig- 
nal strong enough to move beyond the average properties. 
The LBG sample has therefore been split according to op- 
tical magnitud e in a manner similar to that employed by 
E;.aird ct al. (20051) for z ~ 1 galaxies. With the limited 
statistics of only three bins and large error bars we find no 
evidence for a direct correlation between !% and soft X-ray 
flux (Figure|H]and TableEJ, which corresponds to rest-frame 
1800 A and 2-10 keV emission at z ~ 3. This is possibly in 
contrast to that seen for z ~ 1 star forming gala xies where a 
linea r relation between Lx and Luv was found iLaird et alJ 
l2005h . While it is difficult to draw conclusions from this 
data, one interpretation of this result is that at the higher 
SFRs under consideration here, the effects of dust attenua- 
tion are stronger and act to rem ove the direct correlation be - 
tween UV luminosity and SFR llAdelbereer fc Steidel2o'ooll . 
As can be seen in Table |7| the extinction estimates actually 
decrease for the more UV luminous LBGs, although within 
the large errors they are consistent with being constant. This 
suggests that UV attenuation may not be a direct function 
of SFR over a small range in SFRs. 



7 SUMMARY 

Using six ISteidel et al] ^l2003^ LBG survey fields with Chan- 
dra/ ACIS imaging, the X-ray properties of a large sample 
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of UV-selected z ~ 3 LBGs have been examined. The X-ray 
data were used to identify and study AGN within the sam- 
ple, as well as to provide a high-energy perspective on the 
star formation in the galaxies that do not harbour an AGN. 
The main results are: 

(i) 24 LBGs have been detected, approximately doubling 
the number of known X-ray detections. The increased sam- 
ple of X-ray detected LBGs results from ne w analysis of 
Chand ra archive data covering four of the ISteidel et alJ 
l|20o3) survey fields and the doubling of the exposure of the 
Chandra data in the HDF-N, compared to the analysis of 
N02. The AGN fraction in LBG surveys is approximately 3 
per cent, much lower than submm galaxies. 

(ii) Around 1/3 of the X-ray detected LBGs were identi- 
fied as broad line AGN, 1/3 as narrow line AGN and 1/3 as 
normal star forming galaxies at 2 ~ 3. The X-ray luminosi- 
ties of the LBGs ranged from 1.5 x lO"^ to 5 x 10** erg s"^ 
(2-10 keV), therefore spanning Seyfert to quasar luminosi- 
ties. The range in luminosity and the breakdown of spec- 
tral types is similar to that seen in X-ray surveys at lower 
redshift. Deep and comprehensive optical spectroscopy is 
therefore a reasonab ly good way of identifying AGN at z~3 
JSteidel et alJl2002ll . in this case missing about 30 per cent 
of the population. 

(iii) The X-ray to optical flux ratio of the LBGs covers 
two orders of magnitude, ranging from classic AGN to low 
luminosity AGN and starburst values. Two LBGs (HDF- 
M9 and HDF-M35) have low X-ray to optical flux ratio, soft 
X-ray emission, low X-ray luminosity and may be powered 
in the X-ray not by an AGN but by starburst activity. If 
powered solely by star formation, these LBGs have X-ray 
derived SFRs in the range 300 to 500 Mq yr~^ and would 
represent the upper tail of the LBG SFR distribution. Strong 
MIPS 24/im detections and similar FIR derived SFRs lend 
support to the star formation hypothesis. 

(iv) The X-ray spectra were analyzed for each of the 
LBGs with greater than 10 counts in the extracted spec- 
trum. The mean photon index of the X-ray detected LBGs, 
allowing for intrinsic absorption in each of the sources, was 
found to be r = 1.96to:22- 

(v) Significant obscuration was detected in 40 per cent 
of the LBGs. All but one of the narrow line AGN have 
Nh > 10^^ cm~^. Four luminous, obscured type 2 QSO can- 
didates are found after the X-ray luminosities have been cor- 
rected for the observed column densities, although definitive 
confirmation of the type II nature of these objects is impos- 
sible without rest-frame optical Balmer line spectroscopy. 

(vi) A naive conversion from fiux gives a mean luminos- 
ity for the narrow line AGN LBGs approximately half that 
of the QSOs. However, as a group the narrow line AGN 
were found to have significantly higher Nh than the QSO 
group and after correcting the luminosities for absorption, 
the mean luminosity of the AGN and QSO classes were simi- 
lar. The LBGs classified optically as normal galaxies are ap- 
proximately an order of magnitude less luminous than the 
optically identified AGN. The optical light in these weak 
AGN is probably dominated by the host galaxy. 

(vii) Stacking the soft band fiux of 277 undetected LBGs 
in the 2 Ms HDF-N produces a highly significant detection. 
The resulting average luminosity of L2-10 kcv = 2.12±0.39x 
10*^ erg s"^ and X-ray derived SFR of 42±8 M© yr"^ are in 



excellent agreement with that found by N02 using the 1 Ms 
data. The implied UV extinction correction is a factor 4.1 ± 
0.8, consistent with that found in several previous studies. 
Stacking the hard band flux did not produce a detection. 

(viii) Stacking the soft band flux of 226 undetected LBGs 
in the 200 ks GWS produced a marginally signiflcant, 3.3cr 
detection. The derived average flux, luminosity and SFR 
were larger than, but consistent with those found for the 
HDF-N. There is some evidence that the elevated flux is 
due to the presence of AGN within the stacking sample. 

(ix) Splitting the LBG sample into three subsets based 
on observed frame magnitude and stacking the X-ray flux 
produced marginally significant results in each bin, allowing 
the correlation between X-ray and UV emission of LBGs to 
be examined. Unlike for UV-selected galaxies at 2: ~ 1, there 
is no evidence for a correlation between rest-frame X-ray and 
UV fiux. 
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Table 10. X-ray detected LBGs in the lSteidel et alJ i2()0^ fields. Col.(l): Chandra designation (J2000). Col.(2): LBG survey name from lSteidel et al.l l200?t) . except HDF-MD39 which 
is from new LBG survey. CoI.(3): ^ magnitude. Col. (4): Spectroscopic redshift. Col. (5): Positional ofli'set in arcseconds after applying astrometric shift as described in §2.3. Col. (6): Full 
band background subtracted photons detected in 90 per cent EEF area. Col. (7): Lowest false detection probability found for the four detection bands. Probability of 10~* is assigned 
if Poisson probability is less than this value. Col. (8): Full band flux (all fluxes 10~^^ erg cm"^ s"^). Col. (9); Soft band flux. Col. (10): Hard band flux. Col. (11): Rest-frame 2-10 keV 
luminosity in units of 10*^ erg s"^, converted from soft band flux. Col. (12): Hardness ratio, HR = (H — S)/(H -|- S), where H and S are the net hard and soft band counts, corrected to 
on-axis values. Col. (13): Optical spectral classification. Col. (14): Previous identification as X-ray detected 2~3 LBGs: (a) N02, (b) [Basu-Zvch & Scharl 2004, (c) Nandra et al. 2 0oH 
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